Constraints on primordial non- Gaussianity from WMAP7 
and Luminous Red Galaxies power spectrum and forecast for future surveys 
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We place new constraints on the primordial local non- Gaussianity parameter /jvi, using recent 
Cosmic Microwave Background anisotropy and galaxy clustering data. We model the galaxy power 
spectrum according to the halo model, accounting for a scale dependent bias correction proportional 
to fNh/k^. We first constrain /jvl in a full 13 parameters analysis that includes 5 parameters of 
the halo model and 7 cosmological parameters. Using the WMAP7 CMB data and the SDSS DR4 
galaxy power spectrum, we find /jvz, ~ 17ll5^3g at 68% C.L. and —69 < /jvi < +492 at 95% 
C.L.. We discuss the degeneracies between /jvl and other cosmological parameters. Including 
SN-Ia data and priors on Hq from Hubble Space Telescope observations we find a stronger bound: 
—35 < /jvL < +479 at 95%. We also fit the more recent SDSS DR7 halo power spectrum data 
finding, for a ACDM+/ivL model, fNL = -93 ± 128 at 68% C.L. and -327 < /jvi < +177 at 95% 
C.L.. We finally forecast the constraints on /jvl from future surveys as EUCLID and from CMB 
missions as Planck showing that their combined analysis could detect /jvl 5. 



I. INTRODUCTION 

The standard paradigm of structure formation relies 
on the inflation As shown in di], quantum- 

mechanical fluctuations in the scalar field driving in- 
flation lead to primordial density perturbations re- 
sponsible of the large scale structures we observe to- 
day. Although the simplest assumption is that these 
fluctuations were Gaussian distributed [l^, there are 
several inflationary models [ll|-[3| involving the exis- 
tence of a primordial non-Gaussianity. A detection or 
exclusion of non-Gaussianities would hence be of fun- 
damental interest for the understanding the physics 
of the primordial Universe. Several cosmological ob- 
servables and methods can be used to constrain non- 
gaussianities. 

Cosmic Microwave Background (CMB) anisotropics 
provide the most direct method for the detection of 
primordial non-Gaussianity (see e.g. [13) through, 
for example, measurements of the three-point corre- 
lation function [l^ (or equivalently the bispectrum) 
which is non-zero in presence of non-Gaussianities. 
The large scale structure of the Universe is also af- 
fected by non-Gaussianities that may be detected by 
looking at the bispectrum or the trispectrum of galaxy 
distribution. The abundance of galaxy clusters, that 
depends on the tails of the density probability distri- 
bution, is also sensitive to any deviation from gaus- 
sianity. Non-Gaussianity has a direct impact on the 
clustering of dark matter halos by changing their mass 
and correlation function f[l9|-(24|). 
A common way to parameterize primordial non- 
Gaussianities is to introduce a quadratic correction 
to the potential [Hi [111: 



$ = + /jvl(0' -(</'")) 



(1) 



where $ is the primordial potential and is a gaus- 
sian random field. In this case the non-Gaussianity 



is a local type correction whose amplitude is given by 
/jvl- The most recent constraint on /jvl from CMB 
gives -10 < fNL < +74 at 95% C.L. from bispec- 
trum analysis of WMAP-7 years data [23|, improving 
the WMAP5 constrain (-9 < fwL < +111) [II- The 
authors of [2^ used a different estimator applied to 
WMAP5 data, finding /jvl = -13 ± 62 at 68% C.L., 
while [30| , using the needlet bispectrum applied to the 
same data found Jnl = +84 ± 40 at 68% C.L.. 
In |3l|-[3g-[2i|-[33|-[33 it has been shown that a 
quadratic correction to the potential like that of equa- 
tion ([T]) produces a scale dependence in the bias of 
the galaxy clustering with respect to matter distri- 
bution. In particular, a scale dependent term A6(fc) 
arises in the halo bias on larger scales (smaller k) and 
is proportional to /jvl (A6(/c) oc fNL/k"^) hence with 
galaxies being more (less) clustered for pos itive (neg- 
ative) values of /jvl- The authors of [3J| analyzed 
the galaxy power spectrum of luminous red galaxies 
(LRG) of the Sloan Digital Sky Survey (SDSS) to con- 
strain this scale dependence of the bias, putting the 
constraint -21 < /jvl < +209 at 95% C.L.. In the 
same work other large scale datasets have been used to 
constrain non-gaussianity (Quasars, integrated Sachs- 
Wolfe effect data and photometric LRG sample) find- 
ing —29 < /jvl < +70 at 95% C.L. from the combina- 
tion of all datasets. Recently, the authors of [ssj ob- 
tained +25 < /jvl < +117 at 95% C.L. from the com- 
bination of WMAP-7 years data |23|, Baryonic Oscil- 
lations data from SDSS and Two-degree Field Galaxy 
Redshift Survey (2dFGRS) [36| and Supernovae dis- 
tance moduli measurements |37i | with auto correlation 
function measurements of radio sources from NRAO 
VLA Sky Survey claiming a detection of non 

gaussianity at ^ Str. 

In this paper we follow the methodology of [s^l and 
constrain /jvl by looking at the scale dependence of 
the bias in current galaxy surveys data. We imple- 
ment the calculation of galaxy and halo power spec- 
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trum using the halo-model (see [s^ and section|n| and 
we include in it the non-Gaussian scale-dependent cor- 
rection to the bias. We then place constraints on /jvl 
by comparing this model to the Sloan Digital Sky Sur- 
vey (SDSS) [401 galaxy power spectrum data [41| and 
to the halo power spectrum data obtained from the 
luminous red galaxies (LRG) sample We include 
in the analysis the WMAP-7 years cosmic microwave 
background anisotropy data [231 • We also fit the LRG 
galaxy power spectrum data to the same model, in- 
cluding Hubble constant measurements from Hubble 
Space Telescope (HST, [4^) and Supernovae distance 
moduli measurements for the Union dataset [s^l- Re- 
sults are shown in section IHIl We finally forecast the 
power of future Galaxy surveys in constraining non 
gaussianity by generating mock data for galaxy power 
spectrum using specifications of EUCLID [44|| survey 
combined with mock data from Planck [43 satellite 
and showing in section IIVI that the combination of 
data from these experiments could reach the precision 
required to detect even small non-Gaussianities. 



II. HALO-MODEL 



Po{k,z) 
P2h{k,z) = „ , / x 



dMuhaioiM, z)N{M, z)b{M, z)uDM{k, M, z) 



(3) 



where n^aio is the halo mass function |48| . 
UDAiik^ M, z) is the normalized dark matter halo den- 
sity profile in Fourier space, Po{k, z) is the linear dark 
matter power spectrum, b{M, z) the linear bias pa- 
rameter and Ugai is the mean galaxy number per unit 
of comoving volume: 



ngai{z) = j dMnhaio{M,z)N{M,z). 

For low occupied halos {N{M, z) < 1) the exponent 
p of the density profile in the ^ is equal to 1 while 
it is equal to 2 otherwise [s^. To calculate the two 
terms ([2]) and ([3]) it is necessary to assume a form for 
the HOD. We choose the parameterization described 
in [4^ , where the HOD consists of two separated 
contributions for central and for satellite galaxies: 



(Afce„(A/)> - ^Erfc 



\n{M,mn/M) 

cen 



(4) 



In the halo model scenario (see [33 for a detailed 
review) all matter is contained in halos and, as a con- 
sequence, the abundance of halos, their spatial distri- 
bution and their internal density profiles are closely 
connected to the initial dark matter fluctuation fleld. 
Under the assumption that galaxies are formed in 
these halos of dark matter [43 is then possible to 
use the halo model to calculate the statistical prop- 
erties of distribution of galaxies. To this aim, the ba- 
sic quantity is the halo occupation distribution (HOD, 
see [47|) that encodes the information on how galax- 
ies populate dark matter halos as a function of halo- 
mass. The statistical information is contained in the 
two-point correlation function of galaxies or equiva- 
lently its Fourier transform, the galaxy power spec- 
trum. It is hence important to assess the number of 
pairs of galaxies in an individual halo and the number 
of pairs of galaxies in separate halos. The former can 
be shown to be related to the variance of the HOD, 
CT^ (M, z) = {Ng{Ng — 1)) while the latter is the square 
of the mean halo occupation number N{AI, z) = (Ng). 
The galaxy power spectrum is then the sum of the 1- 
halo term describing pairs of objects in the same halo 
and of a 2-halo term for objects in different halos: 
P{k, z) = Pih{k, z) + P2h{k, z). The two terms can be 
written as: 



Pih{k,z) 



gat.-, ^-2) 

dMuhaioiM, z)\uDM{k, M, z)\Pa^iM, z) 



M - ^Mrnir. 

Ml 



(5) 



where Mmim Mi a cen, 7 and a are free parameters of 
the model. In this description the mean occupation 
number of central galaxies is modeled as a smoothed 
step function above the minimum mass M,„i„, while 
satellite galaxies follow a Poisson distribution with 
a mean given by a power low and a cut-off at mul- 
tiple 7 of the minimum mass. This 5-parameters 
model showed a good agreement with hydrodynami- 
cal and N-body simulations and semi-analytic models 

m, [n, M- 

For the halo density profile UDM{k, M, z) we choose 
the shape of the Navarro, Frenk & White profile 
(NFW) [13. The variance of the HOD can be cal- 
culated as in 15511: 



cr(M,z) = N{M,z), 
a{M,z) = /3{MfN{M,z), 



N{M, z) > 1 (6) 
N{M,z)<l (7) 



with l3{M,z) = logio(M/M™„)/Zo5io(Mo/M™™) 
and Mq is the mass at which the mean occupa- 
tion number is equal to 1. This parameterization of 
(t(M, z) has been shown to have a good agreement 
with both semi-analytic models and hydrodinamical 
simulations [5^, (57| . 

The halo mass function is given by the Press & 
Schechter relation ISSl: 



A'PnhaioiM,z) dM 
p M 



(8) 
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FIG. 1: Scale dependent correction to the halo bias ac- 
cording to the (|12p for different values of /jvz,. 



FIG. 2: Best fit galaxy power spectrum calculated with 
the values of table |I] for the fit to WMAP7 + LRG data 
compared with LRG galaxy power spectrum. 



where p is the background comoving density and v 
is defined as the ratio between the critical density 
required for spherical collapse at redshift z {Ssc{z)) 
and the variance of the initial density fluctuation field 
cro{M): V = 5'j.Jz)/i7'^{M). Here we choose the Sheth- 
Tormen model |59j for the shape of vfiv) : 



vf{v)^A{p){l 



1/2 



exp 



-f ) '»> 



with p - 0.3, A[p) - 0.3222 and q ~ 0.75. The linear 
bias h{M, z) is then given by [s^, (60| : 



6(M, z) 



1 



-I 2p/5sc{z) 
5sc{z) 1 + {qvy 



(10) 



A. Non-Gaussian corrections 

The halo model described so far allows the calcula- 
tion of the galaxy power spectrum starting from the 
assumption of gaussian primordial fluctuations. The 
existence of deviations from gaussianity determines a 
correlation between small-scale and large scale pertur- 
bations because of the quadratic correction /nl'^''^ (in 
the case of local non-gaussianities we are considering 
here) that appear in the potential [6l| , [2^ . As shown 
in [sil-fU-H^I the effect of non-Gaussian fluctuations 
on the galaxy power spectrum appears on large scales 
through a scale dependent correction of the halo bias. 
Following we write this correction as: 



Ab{M,z,k) 



3n, 



c^k^T{k)G{z 
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that reduces to: 



Ab{M, z, k) 



c^k^T{k)G{z 



■fNL{b~r)6sc (12) 
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TABLE I: Best fit values and 68% C.L. errors on the 
parameters of our model for W MAP7 -f LRG and 
WMAP7 + LRG + SNe + HST data. The combination 
with SNe ans HST data improves only slightly the con- 
straints. 



where G{z) is the linear growth factor, T{k) is the 
transfer function and the parameter r is 1 if the ob- 
jects equally populate all the halos, that is a good as- 
sumption for the LRG galaxies we use in this analysis 
(see also jsj) or ~ 1.6 for objects populating only re- 
cently merged halos. The information on /at^ is then 
expected to come from the low-fc part of galaxy power 
spectrum because of the k~^ term of the ([12]) {T{k) is 
constant at low wave vectors) . For an even quite large 
value of JnLi for example Jnl = +100, the correction 
on the halo bias is smaller than a 10% for wave vectors 
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FIG. 3: 68% and 95% C.L. contour plots and likelihoods for /nl and other parameters of the model for the fit to SDSS 
DR4 galaxy power spectrum combined with WMAP7, SNe and HST data. The plot shows the degeneracies of Jnl 
with cosmological and halo model parameters. As expected /jvi results to be correlated with matter density and Hubble 
parameter. Strong degeneracies involve also Mmin, ^m, eg and Ho weakening the constraints on these parameters. 



k > (see Fig. [J). The amplitude of the 

correction is proportional to /jvl but also to Hq and 
flrn ■ One should expect hence important degeneracies 
among these parameters that will affect the strength 
of constraints on f^^ ■ 

III. ANALYSIS AND RESULTS 

A. Constraints from Red Luminous Galaxies 
power spectra 

We implemented the calculation of the theoretical 
galaxy power spectrum through the halo model de- 
scribed above and performed a Monte Carlo Markov 
Chain analysis using Cosmic Microwave Background 
data from WMAP-7 years of observations and the 
most recent LRG galaxy power spectrum data (4]| 
available from Sloan digital Sky Survey at a mean 
redshift z ~ 0.35. We fitted these data assuming flat- 
ness of the Universe over a 13 parameters model that 



consists of 7 standard cosmological parameters (the 
physical baryon and cold dark matter densities, the 
ratio of sound horizon to the angular diameter dis- 
tance at decoupling, , the optical depth to reioniza- 
tion, the scalar spectral index, the overall normaliza- 
tion of the spectrum at k = 0.002hMpc~^ and the am- 
plitude of Sunyaev-Zel'dovich spectrum: r2f,/i^, 0.ch? , 
d, r, rig, logiQ 10^'^ As, Asz ) and of the 5 parame- 
ters of the halo-model plus the non-gaussianity pa- 
rameter (logio^min, logigMi, fNL)- In 
what follows we will express the masses in units of 
solar masses. The Markov Chain analysis has been 
performed using the publicly available cosmological 
code cosmoMC jlll suitably modified to include the 
calculation of the halo model and to fit over the pa- 
rameters of the halo model and fNL- The conver- 
gence diagnostic of this code is based on the Gelman 
and Rubin statistic [g^ (also known as i? — 1 statis- 
tic, where R is defined as the ratio between the vari- 
ance of chain means and the mean of variances). The 
results of our fit are shown in table |T] and Fig. [31 
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FIG. 4: 68% and 95% C.L. contour plots and likelihoods for some parameters of our model for the fit to SDSS-DR7 data. 
As noted before /jvi results to be correlated mainly with matter density and Hubble parameter. 



For this model we found weak constraints on the non- 
gaussianity, Jnl ~ i7l^\lg at 68% C.L. and a range 
-69 < /jvL < +492 at 95% C.L. from the combina- 
tion WM API and LRG galaxy power spectrum. 
The best fit power spectrum computed for the values 
of table H] is shown in Fig. [5J as one can see there is 
a slight preference for a non-zero value of /jvl at Itr 
but is largely consistent with gaussian initial condi- 
tions when we consider 2a constraints. These limits 
are weaker than those obtained with a similar dataset 
in (3^ where small scale non linearities are modeled 
with a two parameter k dependent correction. The 
difference is that in our case the uncertainty on /jvl is 
heavily affected by degeneracies with Q.m, Hq and as, 



parameters which are themselves degenerate with the 
parameters of the halo model. We are in fact requiring 
the information on l to come only from SDSS data 
since we are using WMAP data only to constrain cos- 
mological parameters. The LRG data range only for 
scales between 0.01 < k < 0.2hMpc^^ and, as shown 
in Figure [21 on these scales the effect of an even large 
non-gaussianity is small and can be easily confused 
with the effect of cosmological or halo-model param- 
eters. We repeated this fit including both the Hubble 
Space Telescope prior on Hq from |43| and Supernovae 
distance moduli measurements for the Union dataset 
(37! obtaining only a slightly improved constraint on 
Jnl, i-e. —35 < Jnl < +479 at 95%, and other pa- 
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TABLE II: Best fit values and 68% CL. errors on the 
parameters of our model for the fit to SDSS-DR7 halo 
power spectra. 



FIG. 5: Best fit model (solid line) to halo power spectra 
data from SDSS-DR7. We show for comparison the same 
model but with /nl = 0. A negative value of /nl allows 
to have a better fit to the data points in the range (0.03 < 
k < Om)hMpc-^ 



rameters. In Fig [3] we plot constraints on fj^L and 
on the parameters most involved in degeneracies with 
Jnl for the fit to WMAP7 + LEG + SNe + HST 
data. As for the model parameters, we found gener- 
ally higher values for log^Q Mmim 7, and Ucen than 
[49|, but with greater uncertainty, while a results to 
be in good agreement and logj^g Mi has a very large 
uncertainty. These differences may arise because of 
the different dataset and modeling (in [4^ they fit the 
projected correlation function) and parameter space. 
The same figure confirms the expected degeneracy be- 
tween Vim and logioMmin as found also in j49| . due to 
the correlation between halo masses and number of 
galaxies in massive halos (gj, [IHl . The consequence 
of this degeneracy is that our constraint on the matter 
density is Vim = 0.264 ± 0.022 from WMAP7 + LEG 
and hence only slightly better than the constrain from 
WMAP7 alone Vim = 0.266 ± 0.029. 



B. Constraints from halo power spectra 

In this section we constrain /^vl using recent data 
of power spectrum for the reconstructed halo density 
field derived from a sample of LRGs [42| in the seventh 
data release of the SDSS (DR7). The halo power spec- 
trum is more directly connected to dark matter den- 
sity field for a wider k range and this allows to use data 
points in the power spectrum up to A: ~ 0.2/iA/pc~^. 
The main difference with respect to the analysis of the 
previous subsection is that to model the halo power 
spectrum it is not necessary to model the halo occupa- 
tion distribution of galaxies. The halo power spectrum 



in j43| is modeled as: 

Phalo 

(k) ^ Pda7np{k)rDM. damp 

{k)rhaio{k)Fn{k) (13) 

where Pdamp is a power spectrum that account for 
damping of Baryonic Acoustic Oscillations and is cal- 
culated as: 

Pdarapik) = Pa{k)e-^ Pnw{k) (l - c"^) 

(14) 

with Pq being the linear matter power spectrum and 
Pnw is the matter power spectrum with baryon oscil- 
lations removed calculated as in [6^. The value of a 
is chosen fitting the reconstructed halo density field in 
the mock LRG catalogues [Hlfell. The factor F„(fc) 
is a nuisance term defined as: 

K(fc) (l + ai A) ^) (15) 

where bo is the effective bias of the LRG at the effec- 
tive redshift z^ff = 0.313 and k^ = 0.2hMpc-^. The 
likelihood code for halo power spectra froe SDSS-DR7 
is implemented in cosmoMC and performs a marginal- 
ization over the nuisance parameters 60, ai and 02. 
The terms roM^dampik) and rhaio{k) in ([13]) model 
the connection between the non-linear matter power 
spectrum and the damped linear power spectrum and 
between halo and matter power spectrum and they are 
calibrated against numerical simulations (see section 
3 in [i^l for more details). Here we use a modified ver- 
sion of the modeling described so far introducing the k 
dependent bias correction ([T2|) averaged over masses: 

^ J[b{M, z) + Ab(M, z, k)\nhaUM, z)MdM 
^ '^^ / nhaio{M,z)MdM 

(16) 

We then fitted the data varying /atl together with the 
seven cosmological parameters (Vlbh? , Vlch? , 6, r, Us, 
logiQ lO^'^^s, Asz) and minimizing the chi-square by 
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EUCLID+Planck 



Fiducial value 


a 


fNL = +1 


2.23 


fNL = +5 


2.29 


fNL = +10 


2.39 



TABLE III: la errors (68%C.L.)from the combination of 
mock datasets generated for the specifications of Planck 
experiment and EUCLID survey amd for three different 
fiducial values of /nl- 



varying nuisance parameters ai and 02. Our results 
are shown in table and in Figure 01 We also show 
our best fit halo power spectra in Figure [S] Our fit for 
this dataset shows a preference for a negative value 
Jnl = -93 ± 128 at 68% C.L. that allows to have a 
better fit to five point in the observed power spectra 
in the range (0.03 < fc < 0.05)hMpc-^. The uncer- 
tainty on this value remains anyway quite large and 
the 95% C.L. range for Jnl is -327 < /nl < +177, 
hence with a very slight improvement with respect to 
the constraints from previous dataset. For the other 
cosmological parameters we find a good agreement 
with results from for WMAP7 combined with halo 
power spectra of LRG sample. 

We note that, according to degeneracies showed in 
Figure 31 allowing for a possible non-gaussianity im- 
plies an increase of uncertainty on some cosmological 
parameters, namely h and erg with an increase of a 
~ 10% on the la error with respect to KCDM case 
for WMAP7+LRG and an increase of a - 14% on the 
error for ilrh^. 



IV. FORECAST FOR FUTURE SURVEYS 

In this section we consider constraints on /jv l from 
future data generating mock datasets for both CMB 
anisotropy and galaxy power spectra. For a galaxy 
surveys the error on the matter power spectrum can 
be calculated as [631 [69|: 



(17) 



where the effective volume of the survey is given by: 

2 



Veff = V 



nP 



nP+1 



(18) 



and Ak is the width of /c-bins. Here we use specifica- 
tion for a typical future galaxy survey like EUCLID 
[4^ with galaxy number density n ~ 1.6 • 10~^, red- 
shift range < z < 2 and fsky — 0.5. The minimum 
k of the mock dataset is choose to be greater than 
27r/yi/3, while the maximum fc we use is 0.02/iMpc . 
For CMB anisotropy power spectrum we use specifi- 
cation for Planck experiment (45| assuming the noise 



of the 143GHz channel. We explore a A CDM+Jtvl 
model and we choose as fiducial model the WMAP-7 
years best fit for A CDM parameters (23|. For the 
non-Gaussianity parameter we choose three fiducial 
models, /nl ~ +l,+5,+10. Remember that in our 
approach we are using only the information of large 
scale galaxy clustering to constrain non-gaussianity, 
while we use CMB measurements only to constrain 
other cosmological parameters and hence to break de- 
generacies. 

Results for our forecast on f^L are shown in ta- 
ble mil As one can see the combination of accu- 
rate galaxy power spectrum measurements, attainable 
with a survey like EUCLID and Planck CMB measure- 
ments could reach the sensitivity required to detect 
an even small non-Gaussianity, such as /jvl = +5 or 
Jnl = +10, with a confidence level of at least 95%. 
We note that this results is in agreement with other 
forecast for future galaxy surveys (see [T^I for exam- 
ple). Very small non-Gaussianities {Jnl = +1) seem 
instead rather difficult to detect, mainly due to degen- 
eracies with other parameters. Nevertheless in [71| it 
has also been shown that in more complicated mod- 
els (allowing variation of neutrino mass, running of 
spectral index, dark energy equation of state and rel- 
ativistic degrees of freedom) constraints on frfL may 
deteriorate up to ^ 80%. 



V. SYSTEMATICS 

Before concluding we discuss the possible system- 
atics introduced by the assumptions we made or, 
more generally, by the theoretical uncertainties of the 
model. 

First, we have seen that the value of r that appear 
in the ((T2l) may have a value in the range 1 — 1.6. 
We have assumed r — 1 since we are using Luminous 
Red Galaxies that are old galaxies at the center of 
halos. This is a common assumption for this kind of 
analysis (see also [l^l)- Anyway we find that even as- 
suming r = 1.6 the differences in the power spectrum 
with respect to the case r = 1 are very small. Using 
r = 1.6 we find only a slight variation in the with 
respect to r = 1 for the best fit model: Ax^ ^ 0.3. 
The reason for this is that actual data constrain scales 
k > 0.02/iAfpc"^ where the exact value of r is less rel- 
evant. Nevertheless for future low-fc data it may be 
necessary a more precise modeling of the ([T^ . 
A second assumption we made is that the density pro- 
file of the halos is described by the NFW profile. Al- 
though the exact shape of the profiles in the halos is 
still uncertain, this profile has been tested against sev- 
eral numerical simulations and it turned out to be a 
good approximation (39| . Moreover, again, the infor- 
mation on /jvi is coming from k < 0.1hMpc~^ where 
the density profile is constant in Fourier space. 
An important point is the comparison between the 
results from the LRG power spectrum of f4l| and the 
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halo power spectrum from (42|. The first dataset pro- 
vides a galaxy power spectrum while the second a 
halo density field that doesn't require any assump- 
tion about the halo occupation number. In [42| there 
is a large discussion on the differences between these 
two dataset and we refer the reader to this work for 
a complete discussion. Here we remark that the main 
differences are due to the heavy Finger of Gods com- 
pression algorithm used in |4l| to obtain the matter 
power spectrum. This process may cause transfer of 
power from a scale to another, causing consistent devi- 
ations (up to ~ 40% on fc = 0.2hMpc~^) between the 
reconstructed halo density field and the matter power 
spectrum (4^ . For the DR7 halo power spectrum in- 
stead, the halo density field is reconstructed before the 
computation of the power spectrum (see section 2.2 in 
[i^ ) and the deviations between the two are smaller 
than 4%. Also the modeling of the theoretical halo 
power spectrum and of the galaxy power spectrum is 
different. The model of [42| is calibrated on N-body 
simulations and mock data-sets made especially for 
this LRG sample. Moreover the authors of [i^l im- 
posed priors on the nuisance parameters of the model 
based on N-body simulations. For the ga laxy power 
spectrum of [4l| it was used the Q- model 72 1 for the 
non linear part of the power spectrum, marginalizing 
over Q with weak priors. All these differences nec- 
essarily reflects on the cosmological parameters esti- 
mation, including Jnl- The comparison between re- 
sults of the two datasets is made in section 6.1 of (42| . 
Significant differences are found on some cosmological 
parameters from the two LRG datasets only (i.e. not 
including CMB) of the two releases: in particular the 
^Imh values (that enters also in the ((T^)) differ of al- 
most 2a between the two surveys. 
In our work, to fit the DR7 data we are only introduc- 
ing the bias scale dependent correction to the model 
of [i^l, in order to be as much as possible consistent 
with the data compression algorithm of this LRG cat- 
alogue. We ascribe the differences between the results 
of section UlI Al and IIII Bl to the significant differences 
of the data compression process, as noted also in (42| . 
A last issue concern the modeling of power spectra on 
non-linear scales. For DR7 data, the authors of (4^ 
normalize the final halo power spectrum using mock 
catalogues to account for the small offset between the 
N-body and HALOFIT results. 

Concerning the galaxy power spectrum we used to 
fit data from [41|, the P{k) we are using in rela- 
tion ([3]) is the linear matter power spectrum, which 
is well known. The galaxy power spectrum is calcu- 
lated through the halo-model itself. The 5-parameters 
model we are using showed a good agreement with 
hydrodynamical and N-body simulations and semi- 
analytic models [5l|, [Hi, [111, [s^l- Moreover, 
we are marginalizing over the 5 free parameters that 



account for the uncertainties of the model, so that our 
analysis is rather conservative. 

In the forecast section we model the galaxy power 
spectrum relying on the same assumption made above 
(r = 1 for galaxies and NFW profiles) and using the 
same HOD modeling. Our results show the poten- 
tial of a survey like Euclid to detect even small non 
gaussianities and are in good agreement with forecast 
done for the same survey and for a similar modeling 
of the scale dependent bias [t^. It is clear, however 
that the analysis of real data from these future sur- 
veys will probably require a more accurate modeling of 
the galaxy power spectrum and of the scale dependent 
correction in order to not bias the estimated value of 
the cosmological parameters. 



VI. CONCLUSIONS 

We place new constraints on the local type non- 
gaussianity parameter /jv l by looking at the scale de- 
pendence of the halo bias (at small wave vectors) in 
the recent galaxy and halo power spectra measure- 
ments from LRG sample of the Sloan Digital Sky Sur- 
vey. We fit 2006 SDSS power spectra data with an 
halo model consisting of 5 parameters plus 7 cosmo- 
logical paramers and /jv l ■ Our large parameter space 
and the the restriction of the dataset to relatively 
small scales (fc > 0.01hMpc~^) leads to a weak con- 
straint: -69 < /jvL < +492 at 95% C.L.. We show 
and discuss degeneracies with halo model parameters. 
When including both Type la Supernovae and HST 
data the la error on /jvl is reduced of about ^ 10%. 
We use also 2009 halo power spectra data obtained 
from SDSS LRG sample finding a slightly better con- 
straint -327 < /nl < +177 at 95% C.L., again lim- 
ited by the fact that the dataset does not extend below 
fc ^ Q.02hMpc^^. We also forecast the constraints ob- 
tainable from datasets of a survey like EUCLID when 
combined with Planck CMB data, finding that these 
surveys could reach the accuracy required to detect 
even small non-gaussianities as /tvl = +5, thus con- 
firming the power of this method. Finally we discuss 
the possible systematics and theoretical uncertainties 
that may affect the results. 
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